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Abstract 

It has been argued that despite remarkable success, existing random matrix 
theories are not adequate to describe disordered conductors in the metallic regime, 
due to the presence of certain two-body interactions in the effective Hamiltonian 
for the eigenvalues, in addition to the standard logarithmic interaction that arises 
entirely from symmetry considerations. We present a new method that allows exact 
solution of random matrix models with such additional two-body interactions. This 
should broaden the scope of random matrix models in general. 
PACS Nos. 05.40, 72.10.Bg, 05.45. 
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From a phenomenological point of view, random matrix models have proved very use- 
ful in our understanding of a wide variety of physical systems including complex nuclei 
|I| , disordered metals |2| and chaotic systems || . Although the physical systems are very 
diverse, the local statistical properties of the characteristic levels of these systems in the 
bulk of the spectrum turn out to be universal, similar to the well-known universal proper- 
ties of the distribution of eigenvalues of random matrices as proposed originally by Wigner 
. Recently the models have been generalized to include transitions in spectral statistics 
characteristic of metal-insulator or chaotic-regular transitions in finite systems. This 
has opened up the possibility of describing such transitions in this powerful mathematical 
framework, which allows exact evaluation of correlation functions. However, on one hand 
the statistical properties of numerically solved microscopic models with random disor- 
der describing mesoscopic conductors show remarkable agreement with predictions of the 
generalized random matrix theory over a wide range of disorder |B[ ; on the other hand 
there are indications that the appropriate random matrix model for disordered conductor 
is, while highly accurate, not exact in the metallic regime |7], ||. An exact solution || 
for the Fokker-Planck equation describing the probability distribution of the transmission 
coefficients [|1(| shows that the resulting matrix model should include a small correction 



term which apparently destroys the solvability of the model. This correction is responsi- 
ble for a very small correction to the magnitude of the universal conductance fluctuation, 
but at the same time this also resolves a small discrepancy between the random matrix 
result and the perturbative result from microscopic theory JTT]]. While it is not clear how 



important the correction term is to the question of e.g. transition from metal to insulator, 
the fact that the correction exactly reproduces the result of the microscopic theory makes 
it qualitatively non-negligible. It is therefore believed that despite remarkable success, the 



usefulness of the phenomenological random matrix approach for the problem of disordered 
conductors will be severely restricted, if such corrections can not be accomodated within 
a solvable framework. 

In this work we present a new method to accomodate certain type of corrections to 
the standard random matrix models. These corrections are similar to those arising in 
the problem of disordered conductors. The method generalizes the approach based on 
the theory of orthogonal polynomials, and allows exact solutions for physically relevant 
models in terms of known functions. 

The basic ansatz of the random matrix theory is that for a physical system described by 
an NxN matrix X with eigenvalues x n , n = 1, ...N, the joint probabilty distribution (jpd) 
for the ensemble of all random X matrices consistent with given symmetries (hermiticity, 
time reversal etc.) and subject to some physical constraint (e.g. given average density of 
eigenvalues) can be written quite generally in the form || 

P(x u x N ) = J] \x m -x n \ a Y[e- y ^l (1) 



Here a is a symmetry parameter and is equal to 1, 2 or 4 for orthogonal, unitary and 
symplectic symmetries respectively. For example for disordered conductors, a good ansatz 
H is to use the 2iVx2iV matrix X = \[TT^ + (TT^y 1 - 21], where T is the transfer 



matrix characterizing the conductor and / is the unit matrix. The doubly degenerate real 
eigenvalues x are restricted between and oo by current conservation, and directly gives 
the conductance g = J2 n ^ * s useful to describe the probability distribution in terms 
of an effective "Hamiltonian" H of the eigenvalues defined by P = exp(-aH), where 

H{x n ) - - V ln\x m - x n \ + -V(x n ) (2) 

— cv 

m<n 

The repulsive logarithmic 'interaction' term arises from symmetry considerations alone, 



while the confining 'single particle potential' V(x) is the Lagrange multiplier function 



which takes care of the physical constraint [13| mentioned above, and in general depends 



on various physical parameters. For example V(x) = tx, where t depends on disorder, 



describes the disordered metal quite well || 14 



The solvability of the model has so far relied crucially on the fact that the only 
interacting term in (2) is the logarithmic repulsion which arises entirely from symmetry 
considerations; in other words any relevant physical constraint must give rise to only a 
single particle potential. Given this restriction, the universal distributions for nearest 
neighbor spacing or the so called A 3 statistics in the bulk of the spectrum, which we will 
generically call the Wigner distributions [|J, follow from the above jpd when V(x) is taken 
to be linear or quadratic in x. In these cases the potential is strong enough to overcome 
the logarithmic repulsion and the density of levels scale with the number of levels. When 
V(x) is not strong enough, the universality breaks down; in particular for V(x) — > [ln(x)} 2 
for large x there is a transition from the Wigner distribution to an uncorrelated Poisson 
distribution as a function of a parameter ||. Nevertheless, the model remains exactly 
solvable. 

The first hint that a disordered conductor in the metallic regime in higher than one 
dimension may not be exactly described by a simple logarithmic interaction of the above 
form came from attempts to check detailed predictions of random matrix results numer- 
ically [^]. But the nature of the correction needed came from exact solutions |9[] of the 
Fokker-Planck equation satisfied by the transmission eigenvalues in the metallic regime 
[ m | . The solution showed, when mapped to a random matrix Hamiltonian, that the two 
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particle interaction part has the form 

- - ^ ln\x m — x n \ — - ^2 ln\arcsinh 2 (xll 2 ) — arcsinh 2 (x l J 2 )\ (3) 

For x « 1, this reduces to the standard logarithmic repulsion, but for x ~ 1 the 
additional term makes it non logarithmic. The difference is enough to change the variance 
of conductance from the random matrix result | to the microscopic perturbative result j^. 
It is important to establish how significant this small difference is as far as the qualitative 
statistical properties are concerned. But although the existence of such additional two- 
body terms can be understood as arising from some physical constraint that need not be 
of a single particle form, the question of if or how it affects the known random matrix 
results could not be addressed within the current random matrix framework because any 
such additional two-body interaction destroys the existing criterion for solvability, and 
therefore the usefulness, of the model. 

We will show below that with an additional two-body interaction given by a simplifi- 
cation of eq. (3), it is still possible to solve for the model exactly, using a new method. 
While the models constructed are appropriate for disordered conductors, the solvability 
of such models broadens the scope of random matrix theory in general. 

As a first step towards constructing a model that can be solved exactly, and is close 
to a physical model, we approximate the arcsinh function in (3) by a polynomial Sk(x), 
of degree k. In the metallic regime where the model (3) is valid, this should be a good 
approximation. For simplicity and purpose of illustration we will discuss the case Sk(x) = 
x k in detail. We will indicate at the end how the method can, in principle, be generalized 
to arbitrary polynomials. 

We will therefore consider in detail the model described by an additional two body 
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interaction of the form Inix 1 ^ — xjj), which is equivalent to a jpd of the form 

P(xi, x N )= I[(x m -x n )(x k m -x k n ) Y[e- V ^\ (4) 



m<n n 

where k is a positive integer. Note that for k=l the model reduces to the standard unitary 
random matrix ensemble. An exact solution of this model should allow us to understand 
at least the qualitative effects of the additional two-body corrections. 

For the standard logarithmic interaction part we follow the method of orthogonal 
polynomials |4j] and write the product of the differences I\ m<n (x m — x n ) as a Vandermonde 

determinant whose jth column has elements Xx~ , x^ 1 , j varying from 1 to N. The 

determinant remains invariant if we add some linear combination of the other columns with 

lower powers of the x's; the new jth column has elements Y^_i(xi), lj_i(x 2 ), Yj-i{xN), 

where Yj(x) = Yli=objix l is a polynomial in x, of degree j; the coefficients b will depend 
on the choice of the single particle potential V as we will show later. In a similar way, we 
write the correction term x~l m <n{Xm 

second Vandermonde determinant, whose 

jth column has elements Zj-i(xx), Zj-x(x%), Zj-x{xn) where Zj(x) = Z)f=o c jl x i s now 

a polynomial in x k , of degree j; the coefficients c will be determined from the choice of V. 
We now multiply the ith column of each determinant by exp[—V(xi)/2], and interchange 
rows and columns of the second determinant. Eq. (4) can then be written as the product 
of the two determinants, in the form 

P(xx, xn) = Cjs;detK{xi,Xj) (5) 

where is the normalization constant, and 

1 n N-l 



(Vix^ + Vix,)) ^YiMZfa). (6) 
z J 1=0 



K(xi, Xj) = exp 

The reason for writing the jpd as a determinant is the following. Our ability to evaluate 



the n-point correlation function, defined by 



7V! f f 

Rn(xi, x N ) = _ J ... J dx n+1 ...dx N P N {x u x N ) (7) 

depends on our ability to integrate the jpd over arbitrary number of variables. These 
integrals can be done in a very simple way [|J if the jpd can be expressed as a determinant 
as in eq. (5), provided the following two conditions are satisfied: 

(i) J K(x, x)da(x) = constant 

{ii) [ K(x, y)K(y, z)d/i(y) = K(x, z) (8) 



where d\x is a suitable measure. This is where the restriction of the standard logarithmic 
interaction, equivalent to the case k = 1, comes in. For k = 1, the polynomial Z(x) is the 
same as Y(x), and they can be chosen such that they form an orthonormal set p(x) with 
respect to the measure exp[— V(x)]dx, i.e. 

Pm(x)p n (x)e~ v{x) dx = 5 mn . (9) 

Then the above two conditions in (8) follow from the orthogonality and normalizability 
of the polynomials. The additional two-body interaction forces the polynomials to be 
distinct, destroying the orthogonality. We will now show that even for distinct poly- 
nomials Y(x) and Z(x), the above two conditions in (8) can be satisfied under some 
conditions, making it possible to obtain exact solutions for the correlation functions for 
these generalized models. 

Let us choose the coefficients b and c in such a way that the polynomials Y and Z 
satisfy the following: 



J Y n (x)x kj e- Vix) dx = 0, j = 0, 1, ....n - 1, 



0, j = n, 

J Z n (x)x j e- V{x) dx = 0, j = 0, 1, ....n - 1, (10) 
^ 0, j = n. 



It can then be shown JT5[ that the two polynomials form a "biorthogonal" pair, defined 
by 

J e- y ^Y n (x)Z m (x)dx = gj mn (11) 

We will always choose an overall multiplicative factor such that g n — 1, i.e. the polyno- 
mials are normalized. Clearly the two conditions in (8) are then satisfied again, 



N-l 



1=0 



where we have used the normalization, and 



f K{x,x)dx = f e~ v{x) Yi{x)Zi{x)dx 
J J 1=0 

N-l . 

= E j e^Y^Z^dx = N (12) 



N-l 



K(x,y)K(y,z)dy = / e 4[nWl e 4[WW] £ Y^Z^Y^Z^dy 



3,1=0 

N-l , 

= e-i^W+^Wl £ Y j {x)Z l {z) / e-vMYMZjWdy = K(x, z) (13) 

5.1=0 

where we have used the biorthogonality of the polynomials. Given these properties, the 
integration over N — n variables x n+ \, ....xn in the jpd can be explicitly carried out ||, 
and we obtain 

R n = det[K(x i ,x j )] iij=li _.. n (14) 

where the kernel K(xi, Xj) is given by eq. (6). In particular the density is given by K(x, x) 
and the two-particle kernel from which the nearest neighbor spacing distribution or the 
A 3 statistics can be calculated is given by K(x,y)K(y,x). 
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The model then is exactly solvable if for a given choice of the single particle poten- 
tial V(x) the corresponding biorthogonal polynomials can be obtained. The physically 
interesting model that already gives a very good description of the metallic regime of a 
disordered conductor in the k — 1 limit is given by V(x) = x, < x < oo. The model 
is exactly solvable in terms of Laguerre polynomials. For arbitrary k, the correspond- 



ing biorthogonal polynomials has been explicitly calculated by Konhauser | 16fl . [For the 



simplest non trivial case k = 2, these are the polynomials introduced by Spencer and 



Fano \17\ to study penetration of matter by gamma rays, and studied later by Preiser 



18| .1 Therefore using the new method the exact solution can be immediately written 
down in terms of these Konhauser biorthogonal polynomials. It has been argued that an 
appropriate generalization for all disorder, in the k = 1 limit, is given by the choice [[H 



V(x;q) = ^ln[l + (l-g)g n x], < q < 1. (15) 

n=0 

As q — > 1~, V(x) — > x, and one recovers the metallic regime, while increasing disorder cor- 
responds to decreasing q. This model is exactly solvable in terms of the g-Laguerre poly- 
nomials. For arbitrary k, again the corresponding biorthogonal polynomials are explicitly 



known ||19|| , and the exact solution can be written down in terms of these "q-Konhauser" 
biorthogonal polynomials. The detailed properties of these solutions are under investiga- 
tion. 

It is also possible to consider a more general form of the jpd given by 

P(xi, x N )= l[[r(x m )-r(x n )}[(s(x m )-s(x n ) \{e- v ^ . (16) 

m<n n 

where r(x) is a polynomial of degree h and s(x) is a polynomial of degree k. Defining Y 
and Z as polynomials in r(x) and s(x) respectively, the above method should be applicable 



if conditions (10) are replaced by (15 



/ 



e- vix) Y n (x)[s(x)] j dx = 0, j = 0, 1, ....n - 1, 



^ 0, j = n, 



I 



e- v(x) Z n (x)[r(x)} j dx = 0, j = 0, 1, ....n - 1 



(17) 



^ 0, j = n. 



The case considered before is a special case where r(x) = x and s(x) = x k . Note that 
writing r(x n ) - r(x m ) = (x n - x m )r(x m , x n ) and s(x n ) - s(x m ) = (x n - x m )s(x m , x n ), we 
can write the jpd in the form 



which may allow more physically interesting models to be solved exactly, if the corre- 
sponding biorthogonal polynomials are known. 

In summary, we present a new method to accomodate certain two-body interactions in 
random matrix models, particularly appropriate for the problem of transport in disordered 
conductors. We show that correlation functions can be written down explicitly in terms 
of known biorthogonal polynomials. The approach should broaden the scope of random 
matrix models in general. 

I am grateful to Mourad Ismail for many discussions, and in particular for bringing 
reference [19] to my attention. 



P(xi, xn = Y[(x m -x n ) 2 Y[ r(x m ,x ri )s(x m ,x n ) Y[e 



(18) 
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